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Abstract. A linear response framework is set up for the evaluation of collective 
excitations in a confined vapour of interacting Bose atoms at finite temperature. 
Focusing on the currently relevant case of contact interactions between the atoms, 
the theory is developed within a random phase approximation with exchange. 
This approach is naturally introduced in a two-fluid description by expressing the 
density response of both the condensate and the non-condensate in terms of the 
response of a Hartree-Fock reference gas to the selfconsistent Hartree-Fock potentials. 
Such an approximate account of correlations (i) preserves an interplay between 
the condensate and the non-condensate through off-diagonal components of the 
response, which instead vanish in the Hartree-Fock-Bogolubov approximation; and 
(ii) yields a common resonant structure for the four partial response functions. 
The theory reduces to the temperature-dependent Hartree-Fock-Bogolubov-Popov 
approximation for the fluctuations of the condensate when its coupling with the 
density fluctuations of the non-condensate is neglected. Analytic results are presented 
which are amenable to numerical calculations and to inclusion of damping rates. 



PACS numbers: 03.75.Fi, 05.30.Jp 
1. Introduction 

The achievement of Bose-Einstein condensation in trapped vapours of alkaU atoms 
0, H, H is allowing one to probe properties of an inhomogeneuos dilute Bose gas 
which were not accessible earlier. In particular, the excitation of shape deformation 
modes in a regime for which there is no detectable non-condensate fraction jj, ^ 
has allowed tests of theoretical studies of the collective excitations based on a mean- 
field description of the condensate at zero-temperature An extension of the 
theory to finite temperature has also been given within the Hartree-Fock-Bogolubov 
approximation (HFBA), with an explicit calculation being performed for the case of 
isotropic harmonic confineme nt p3[ . 

Very recently Jin et al. |14[ have extended the experimental study of low-lying 
collective excitations to finite temperature. Their measurements were carried out on 
two modes with different angular momenta and show very striking dependences on 
temperature for both the mode frequencies and the damping rates. These authors 
stress the importance of an interplay between condensate and non-condensate as a 
potentially dissipative process which is not included in the usual mean-field theoretical 
models. 
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In this work we present what may be regarded as the first step in transcending 
a mean-field treatment at finite temperature. This is taken through inclusion of 
correlations via a random phase approximation with exchange (RPAE). Our approach 
is inspired to RPA theories of the linear density response in homogeneous quantal 
, |l6[ |l^, which construct the response of the interacting system to weak 
external perturbations via the response of an ideal reference gas to selfconsistent 
Hartree potentials. However, in the experimentally relevant situation for the confined 
Bose-condensed systems we need to include exchange in the selfconsistent potentials 
set up by density fluctuations in both the condensate and the non-condensate. The 
requirement of consistency between the equilibrium state and the fluctuations around 
it then leads us to choose a Hartree-Fock g reference system. This choice also 

permits us to regain the HFBA limit at finite temperature when the non-condensate 
fluctuations are neglected. 

The layout of the paper is briefly as follows. Section 2 sets out some general 
aspects of linear density response theory for a two-component system representing an 
inhomogeneous Bose system in the presence of a condensate. After an evaluation in 
Section 3 of the density response functions in the HFBA, Section 4 introduces the 
finite-temperature RPAE for the case of contact interactions between the particles 
and discusses a general method of solution of the RPAE integral equations in order 
to display the resonant structure of the partial density response functions. Section 5 
derives the RPAE equations by an equation-of-motion approach and elaborates the 
choice of the reference gas, leaving the comparison with the HFBA theory to an 
appendix. Finally, Section 6 contains a summary of the results and a discussion of 
further developments. 



2. Linear density response functions 

We consider a confined system of Bose particles described by the Hamiltonian Hg, 
which is subjected to an external perturbation represented by a scalar potential 
Up{r,t). The perturbation Hamiltonian is 

H{t)^ J d^rh{r)Up{r,t) (1) 

where n(r) = ^'''(j.^vi^^i-) jg tj^g particle density operator and ^(r) is the field operator. 
The density change to first order in the perturbation is given by 

t r 

'3' 



Sn{r, t)^i I dt' d^r' ([n(r', t'), n(r, t)])eg Up{r' , t') 

dt' f d^r'xnn{r,r';t-t')Up{r',t') , (2) 



defining the linear density response function Xrm(r, r';t — t'). In equation we have 
set fi{r,t) — exp{iHQt) h{r) cxp(—iHot) {fi — 1) and have assumed that the system 
was initially at equilibrium. {0)eq denotes the statistical expectation value of the 
operator O in the equilibrium ensemble. 

In the presence of a condensate, the field operator 'I'(r, t) is split into the sum of 
the condensate wavefunction $(r,t) = (5'(r,t)} and of the field operator ^'(r, i) for 
the non-condensate. The density response function is then given as the sum of four 
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components, 



Xnn{r, r'; t) = Xcc(r, r'; t) + Xcn(r, r'; + Xnc(r, r'; t) + Xnn(r, r'; t) (3) 



where 



+ 16 



Xcc(r, r'; i) = ($:,(r')$e,(r, t)( [*(r'), (r, t)] }eg + ff.c.) 

Xcn{r, r'; i) = ($:^(r, t)( [*t(r')^'(r'), *(r, <)] ) 
X^.e(r, r'; = ($:^(r')( [*(r'), *^(r, t)^'(r, t)] } 



eg 



eg 



H.C , 



and 



(4) 

(5) 
(6) 

(7) 



Here, H.c. denotes the Hermitian conjugate. These equations define a two-by-two 
matrix of response functions for the inhomogeneous two-component fluid. Within the 
linear response regime the condensate wavefunction in is the equilibrium one, 

with $eg(r, t) = $eg(r) exp(— i^i) and fi being the chemical potential. 

A more standard approach in linear response theory for a Bose fluid in the presence 
of a condensate would actually consider applying an external scalar potential Up(r,t) 
coupled to the total density and a gauge-breaking field rip{r,t) coupled to the field 
operator vE'(r) through the perturbation Hamiltonian 



H{t)^ d\ C/p(r,t)#t(r)^(r)+,,p(r,t)#t(i.)+^*(r,i)*(r) 



(8) 



(for the case of a homogeneous fluid see Griffin jl^l). The change in total density is 
in this case given by 

5n(v,t)^<^,,{v,t)5{¥{v,t))+^:^{v,t)5{i!{v,t))+5n{v,t) . (9) 

Explicit evaluation of the first-order changes induced in the field operator and in the 
non-condensate density in equation (^) shows that this formalism yields the same 
results as in equations provided that 

7yp(r,t) = {/p(r,t)$eg(r,i) . (10) 

We shall in the following adopt the two-fiuids framework set out in equations (|^)- 
(0), which will allow us to treat the condensate and the non-condensate on the same 
basis. However, at some later stage it will be convenient to adopt the expedient 
of considering the application of different external scalar potentials to the two 
components of the fiuid, for the purpose of separately determining the four partial 
density response functions in equation (0). 



3. Density response in the Hartree-Fock-Bogolubov approximation 

In the HFBA the field operator for the non-condensate is expanded in terms of 
normal modes, 



(11) 
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where aj and ctj are the quasiparticle annihilation and creation operators, Uj and 
Vj are the quasiparticle amplitudes and Ej are the excited state energies. These 
operators satisfy the canonical commutation relations and their average on the system 
at equilibrium is given by the Bose factor f{Ej), 

{a]ak)eg = Sj^k (e'3(^^-^) - l)"' = fiEj)Sj^k (12) 
with p = [kBTy^. 

The partial density response functions are easily evaluated by using quasiparticle 
amplitudes given by the HFBA, with the following results: 

Xcc(r,r';t) = ^e{t) \^l^[r')^,,{v,t) ^ {u,{r')u*{v)e^^^' - v*{v')v,{v)e-^^^') - H.c. 

+ <i>e, (r')$e, {v,t)Y^ {u* {v')v, (r)e-^^ * - v, {v')u* (r)e^^^ *) - iJ.c. i (13) 



Xfis(r,r';t) = ie{t) <| ^hjk 

3,k 



u* (r')wfe (r') (Wfe (r)uj (r) + Vj {r)vl (r)) 
e^i^. (r')^;* (r') (li* (r)?., (r) + vu {v)v* (r))" 

-.(i^. (r') (r)x;; (r) + u* (r)vt (r)] ) - iJ.c. Wl4) 



and 



Xfic(r,r';t) = 
Xcn(r,r';i) = 



(15) 
(16) 



Here, hjk = fiEj) ~ f{Ek) and /-,a; = 1 + f{Ej) + f{Ek). These expressions are 
fully determined since $eg(r, t) satisfies the Gross-Pitaevskii equation with eigenvalue 
/i, while Uj(r) and Vj{r) at energy Ej satisfy the Bogolubov equations in the Popov 
approximation p^ . 

As is evident from the structure of equations (^3|)-(|l6|), the HFBA treats the 
excitations of the condensate independently of those of the non-condensate. We briefly 
point out below, for the sake of completeness, how results already known for the 
condensate can be recovered from equation (^3|). We then discuss the predictions 
made by the HFBA in equation (p^ for the excitations of the non-condensate. 



3.1. Oscillations of the condensate 

The density change of the condensate is given in the HFBA by the response function 
Xcc and consists of two terms describing the excitation of particles out of the 
condensate (first term in equation (p^)) and their de-excitation (second term in 
equation (p^). 
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With the notation 

{{uj\Up\^,g)) = I f dt' f d\' e'^^*'zi*(r')C/p(r',t')*e,(r') (17) 

J — OO J 

etcetera, we have for the excitation term the expression 

5nr{v,t) = -$:,(r,i)^ K-(r)e-^^*(^,|C/,|<&e,)) -t';(r)e^^^*((t.;|C/,|<J>e,))] . (18) 

In equation ( p^ ) we recognize the typical structure Sn'^!''^ oc $eq("j ~ "^j) Q ^"^^ the 
existence of selection rules such as used by Dodd et al. ||2^ to evaluate the excitations 
of the condensate in the presence of vortices. With the same notation the de-excitation 
part is given by 

Jr.f--(r,i) = -$:^(r,t)^ [i;;(r)e'^^*((ci>e,|L/,|^.,)) -z.,(r)e-^^*(($e,|f/p|«;))] .(19) 

j 

The excitation energies Ej — fi correspond to resonances in Xcc(r, r'; w). In relation 
to experiments ^, ||, |l^ , the HFB A results account for the observed eigenfrequencies 
of harmonically confined condensates at zero temperature |]^- |pl| , but do not seem to 
account for their observed temperature dependence jl^ . Of course, as a mean field 
theory the HFBA docs not include damping of the modes. 



3.2. Oscillations of the non- condensate 

The oscillations of the non-condensate are described by the h — h response function 
given in the HFBA by equation (|lj). It is seen from that equation that these 
oscillations are of two types, corresponding to the excitation of single quasiparticles 
(first and second term in equation (^4|)) and of pairs of quasiparticles (third term in 
equation ([ij)). The density change of the non-condensate associated with the former 
type of excitations is 

Sh^'Hr,t)=Y,ifiE,)-fiE,))[{uUr)u,ir)+v,ir)vUr)) {{u,\U,\u,)) e-'^^^-^^^' 

- («;(r)u,.(r)+t;,(r)t;;(r))((«*|[/,|t;:))e*(^^-^'=)*) . (20) 
For the two-quasiparticle excitation process we have 

5n^^\v, t) ^ ^(1 + f{E,) + f{Eu)) ({ul{v)v*{v) + u*{v)vl{v)) {{v*\UM) 

- {vu{v)u,{v)+v,{v)u,{v)) ((u,|{/p|z;fc))e-(^^+^'=)*) . (21) 

From equation (|2^) we remark that the two-particle process in the HFBA vanishes in 
the dilute-gas limit at zero temperature. 

In the recent experiments of Jin et al. [Q two discrete eigenfrequencies of the non- 
condensate have been observed. For temperatures higher than Tc these excitations are 
to a good approximation described by equation ( pO| ) evaluated in the ideal-gas limit. 
This is perhaps not surprising, since the non-condensate is in a very dilute regime. 
The same approximation seems to hold for the eigenfrequencies of the non-condensate 
even in a limited range of temperatures below Tc, in which the HFBA fails to account 
for the eigenfrequencies of the condensate. 



6 



4. Random phase approximation 

Here and in the following sections we develop an approach invoking a random phase 
approximation with inclusion of exchange (RPAE). Its main feature is that it treats 
the dynamics of condensate and non-condensate fluctuations on the same footing. 

The need to include the dynamics of non-condensate fluctuations into the theory 
comes from the aforementioned experimental observations of Jin et al. of both 
condensate and non-condensate normal modes. In this spirit, Zaremba et al. ^ have 
recently pioneered an approach which develops hydrodynamic equations including a 
time-dependent non-condensate density. By allowing the non-condensate to fluctuate, 
it is possible to satisfy the generalized Kohn theorem |^ and to obtain the analogue 
of "second sound" as an out-of-phase fluctuation of condensate and non-condensate 
densities. 

The RPAE works in the linear response regime and does not include the damping 
mechanism due to nonlinear interactions between the condensate modes (these 
large amplitude effects are instead contained in the time-dependent Gross-Pitaevskii 
equation |^). However, the RPAE keeps into account the coupling between 
the dynamics of the condensate and that of the non-condensate. This mechanism 
is thought to be important for an explanation of the behaviour of the observed 
frequency shifts with temperature. We support this by noticing that these shifts 
become significant when the temperature of the system becomes of the order of the 
chemical potential, that is when the non-condensate cloud starts being significantly 
populated. 

In order to define the RPAE for the specific system of an inhomogeneuos Bose- 
condensed gas, we introduce at this point a matrix of "proper susceptibilities" X) 
by analogy with a well known approach to homogeneous fluids [ p^ [T^ . This 
matrix relates the density changes of the condensate and the non-condensate to the 
selfconsistent Hartree-Fock potentials. We write the latter in the form 

UHp{r,t) = Up{r,t) + a6nc{r,t) + 2aSii{r,t) (22) 

and 

Ufjp{r,t) = U^{r,t) + 2aSnc{r,t) + 2aSn{r,t) . (23) 

We are using different external scalar potentials for the two components of the fluid, 
as a device to determine the four partial density response functions. Furthermore, 
we are considering the special case of contact interactions between the particles, with 
a — 4:Trh^a/m where a is the scattering length and m is the particle mass. 

The RPAE follows naturally by evaluating the diagonal elements of x from a 
suitably chosen "reference" gas and by setting its off-diagonal elements equal to zero. 
As we shall see explicitly below, the latter approximation on x still preserves the 
off-diagonal elements of the response matrix of the fluid. 

After Fourier transform with respect to the time variable the RPAE equations thus 
read 



6nc{r, io)^ J d^r' x°(r, r'; u) ([/^(r', to) + aSndr' , u) + 2a6n{r' , uj)) (24) 
and 

Sfi{r, Lu)= f d^r' x~ (r, r'; uj) (f/^(r', uj) + 2aSnc{r' , lu) + 2a5n{r' , uj)) . (25) 
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In these equations and Xa ^-re the density response functions of the condensate and 
non-condensate components of the reference gas. 

The choice of the reference gas will be discussed in the next section. Here, instead, 
we want to display the general structure of the solution of equations (|2^) and (p5|). 
A formal solution can be obtained by expansion into a complete orthonormal set of 
basis functions, so as to convert the problem into one of matrix inversion, and by 
comparison of the results with equations and (j^) (extended to the case of two 
different scalar potentials). The results for the density response functions in matrix 
notation are as follows: 

Xcc - (1 - 2axl) [(1 - «Xc) (1 - 2«X^) - 4«'x°X°] X° (26) 
Xcn = 2ax" [(1 - «Xc) (1 - 2ax^) - ^a^'lxl] (27) 



2ax^ [(1 - ax°) (l - 2ax°) - ^a\'cXl] X° (28) 



and 



Xnn = (1 - aXc) [(1 - «X^) (1 - 2ax5) - 4a\°xl] xl (29) 

The appearance of a common "denominator" in the expressions (p6|)-(p9|) reflects the 
fact that the RPA is in general leading to collective excitations of the system as a 
whole. 



5. Equation-of- motion approach to the random phase approximation 
with exchange 

At this point we need to specify the reference gas to be adopted in the RPAE for 
the inhomogeneous Bose fluid. The reference system is to be described by a set of 
single-particle orbitals and energy levels. In the case of a homogeneous Bose fluid, 
translational symmetry forces the choice of a plane-wave set of orbitals. Furthermore, 
contact interactions in fc-space merely add a constant term to the energy. As we 
shall elaborate in the discussion below, we believe that an optimal choice of the 
reference system for the inhomogeneous fluid in the RPAE is provided by the Hartree- 
Fock theory (HFA). This choice (i) allows consistency between the treatment of the 
equilibrium state and that of the fluctuations around it, and (ii) yields results which 
in the appropriate limit reduce to the HFBA treatment for the condensate at flnite 
temperature in the Popov approximation (see Appendix A). 

We begin by examining the equations of motion for the condensate wavefunction 
$(r,f) and for the non-condensate field operator ^{r,t) within the Hartree-Fock 
approximation. These are [ p5| : 

t — '^={-— + Ve + a\'i'\^ + 2an]^ + m<S> (30) 
ot \ 2m J ' 

and 

*9t*" ( -^ + K + 2a|$p + 2an J * + . (31) 

In the above equations we have defined the non-condensate density as n = In 
the case of weak external perturbations the equations can be linearized to read 

( + — - V; - a|$egp - 2ano j 5^ = {Up + aSuc + 2aSfi) $eg (32) 
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and 




eq 



(33) 



where we have set Sric — (5|$p. By operator inversion we find 
5$ = G^^ {U^ + aSn, + 2a6n) $eg 



(34) 



and 



eq 



(35) 



With Go we have indicated the Green's functions in the Gross-Pitaevskii (GP) or 
Hartree-Fock (HF) approximations and have left implicit the convolution integral 
required in the inversion. 

The structure of equations (|3^) and (|3^) can now be compared with that of 
equations (|2^) and (^5|). It is evident that, in order to obtain a consistent treatment, 
we need to construct the reference gas susceptibilities Xc ^^'^ Xfi with single particle 
orbitals which are eigenvectors of the Gross-Pitaevskii operator and of the Hartree- 
Fock operator, respectively. Although both the condensate and the non-condensate 
are thereby being consistently treated in the Hartree-Fock theory, this will necessarily 
lead to the use of different single-particle orbitals and energy levels in the evaluation 
of the two reference-gas susceptibilities. This is a peculiarity of the Bose-condensed 
system, which comes from the fact that there is no exchange in the condensate. 

We come now to a more detailed derivation of the RPAE equations, using the 
equation of motion for the density matrix as in the treatment of the Fermi fluid given 
by Singwi et al. |2^. In a Bose-condensed system the density matrix of the non- 
condensate is defined as 



Of course, the factorized form of Pc shows that this is not a true density matrix, but 
merely a convenient device allowing a "two-fluid" visualization. 

In the occupation-number basis |25| we express the non-condensate field operator 

as 



p(r. 



■ r',i) = (^'t(r,i)*(r',i)) - $*(r, i)$(r', t) = (^^(r, t)^(r', i)) . (36) 




(37) 




(38) 



and the condensate wavcfunction as 




(39) 



The functions introduced in ( p6| ) and ( pT] ) become 



(40) 



and 



{i\Pc\j) = z*z. 



(41) 



5.1. RPAE equation for the non- condensate 

The equation of motion for the non-condensate density matrix in the HFA is 
d . 



X! i^^r^Prj - Ptrhrf + U^rPr] " PirU^^j) (42) 



where h^^^ is the HF single particle Hamiltonian expressed in the occupation number 
basis, 

hf,'' = + V,\:i) +2Y,mv\lj){zlzi+~pu) . (43) 

By linearizing the time-Fourier transform of equation (^) in the HF single particle 
basis (eigenfunctions V'j(r) with eigenvalues Cj) we obtain 

[lo - (e, - £,) + ir^)6~p,, = (/(e,) - /(e,)) f f/^ + 2 ^(*fc|«|/j)('^(4^0 + ^Pki) ) (44) 



fc,; 



where we have set pj^ — fi£j)Sij for the equilibrium density matrix, with i and j 
different from zero since it describes the excited states, and we have introduced a 
factor 77 = 0^ to allow for an adiabatic turning on of the external potential. 

In order to come back to the real space description we recall that the non- 
condensate density has been expanded as 

(5n(r,w) = ^i/..(r)V'*(r)5/5y , (45) 

while the condensate density is given by 

6n,{r, ^) = E V'.(r)V'*(r)<5(z;z,) ; (46) 

moreover for contact interactions we find 



Y^{ik\v\lj)iSizlzi)+Spki) = aJ2 / d^'-r'^lVk^jM^izlzi) + 5pki) 

kl •' 

a I d^r'i;*^j{6nc + Sfi) . (47) 



kl kl 



If we now express the reference response function for the non-condensate as 

X°(r,r',c.)= ^^r^'^^'f ■ V';(rOV',(r)V:(r)^fc(rO , (48) 



by using the relations (|4E| - [48D in the sum over all the modes of (^) we obtain 

^ V;(r)^.(r)<5py = / dV'x?(r,r',c.) ([/; + 2a{Sn, + Sfi)) . (49) 



.0 



This is equation (pa), with Xfi given by equation (|48|). 
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5.2. RPAE equation for the condensate 

The derivation of the RPAE equation for the condensate follows the same lines as 
that for the non-condensate, with some differences that we wish to remark. The HFA 
equation of motion for pc in the occupation number basis is 

^'O^Pij ~ 5Z ip-fr' Prj ~ Pir^^f + U'^^p'^j.j — p^^C/^^) (50) 
r 

where /i,^^ is the Gross-Pitaevskii single particle Hamiltonian 

hf,"" - - 1^ + ^e|j) + Y.'^^k\v\lj){zlzi + 2pki) ■ (51) 

k,l 

One must choose the eigenstates of the Gross-Pitaevskii operator ('/)i(r), say, with 
ei denoting the corresponding energy levels) in order to diagonalize the equation of 
motion for the condensate. 

Within the same scheme the equilibrium density matrix for the condensate is 
(z*Zi)o = Si^oSj,o Nq, where A^o is the condensate fraction, depending on temperature 
as determined from a thermodynamic treatment p7| , The linear fluctuation of 
Pc then is Sp'ij = SpQ^Si o + ^P'ioSj,o- With these definitions, the linearized Fourier 
transform of equation ( po|) is: 

{uj - (cj - p) + iv)^PM Sjfi + {uj + {cj - p) + iv)Spoj Sifi = (52) 

= NoiSj,o - S.,,o) I C/^- + Y.{ik\v\lj)iSizlzi) + 2Spki) 

It is now easy to go back to the real-space equation by using the expression 
of the condensate density, 

6nc = ^ {z^(t)Q(j)^ + Z*(t)Q(j)*) (53) 

and of the condensate reference susceptibility, 

X°(r,r',^) = ( ^_(e^.'!°^)+,/ o(rO05(r)0,(r)0;(rO 



3^0 



MrWoir')^3ir'WAr)] ■ (54) 



UJ + ej — p + irj 



In summary, equations (^8|) and ( p4D give the susceptibilities of the reference gas 
in the RPAE. They involve the Hartree-Fock and the Gross-Pitaevskii orbitals and 
energy levels, respectively. It can also be shown that from equations (42) and (^) the 
HF theory leads to hydrodynamic-like equations for the particle densities and current 
densities, of a form which ensures that the generalized Kohn theorem 22 is satisfied 
(Minguzzi and Tosi, to be published). 



6. Summary and future directions 

We have in this work set up an RPAE formalism for the collective dynamics of an 
inhomogeneous Bose gas at finite temperature. The system is treated within the linear 



11 



response theory as a "two-component fluid" , the components being the condensate and 
the non-condensate. 

The RPAE assumes that the fluid responds to the self-consistent Hartree-Fock 
potentials via the response functions of a Hartree-Fock gas. This ensures consistency 
between the treatment of the equilibrium state and that of its deformations induced 
by external perturbations . It also ensures that the HFB-Popov theory is regained 
when the coupling with the fluctuations of the non-condensate is neglected. 

In essence, what we have proposed is a practical way to go beyond the HFBA 
by including the coupling with non-condensate fluctuations into the description of 
the Bose-condensed system. A microscopic derivation of the RPAE equations by 
an equation-of-motion approach has helped to bring more physical insight into the 
approximations that we have proposed. As we have pointed out, the RPAE gives 
some account of correlations between condensate and non-condensate by preserving 
non-zero off-diagonal elements in the susceptibility matrix. 

The strategy one would adopt to use the RPAE in actual calculations may be 
briefly mentioned. One has to (i) determine the single particle orbitals and energy 
levels of the reference system as functions of temperature, in order to build the 
susceptibilities given in equations ( ^ ) and (54); (ii) solve the integral equations ( p^ 
and (|5|), thus obtaining the four partial response functions by comparison with the 
definitions (^)-(0); and (iii) search for the resonances of the response in the complex 
frequency plane. These give the excitation frequencies of the system as a whole. 

Some additional simplifications seem possible in treating the experimentally 
relevant case of dilute Bose-condensed gases: (i) the role of the interactions in the 
reference gas of the non-condensate is presumably minor; and (ii) in view of the 
resonant structure of the denominators in equations (^) and (|5^) , the search for low- 
lying collective modes should mainly involve the single-particle levels which lie in the 
same frequency range. 

Let us now comment on the insights on further progress that can be drawn from 
the development presented in this work. The issue of the mechanisms of damping 
of the collective excitations is of main interest here. It is relevant to recall that 
damping of sound waves in a homogeneous Bose fluid was obtained in an RPA 
treatment by Szepfalusy and Kondor fl^ . However, in that case the condensate 
is coupled to a continuum of levels. Instead, the single particle orbitals for the gas 
in harmonic confinement have discrete energy eigenvalues, with spacings which are 
not much smaller than the low-lying collective excitation frequencies. Broadening of 
the single-particle levels can be understood as due to the interplay between thermal 
fluctuations and inter-particle interactions and could be very simply included in the 
RPAE equations in a phenomenological fashion by replacing the parameter rj in 
equations ( ^8| ) and ( ^4[ ) by temperature-dependent relaxation rates. Such an approach 
was proposed in early work of Mermin to take collisions between particles into 
account in the RPA dielectric function of the electron gas. 

At a more fundamental level, important damping mechanisms for the collective 
excitations should be associated with those couplings between the condensate and the 
non-condensate which in a linear response framework are described by dynamic off- 
diagonal proper susceptibilities. These could be taken into account in a perturbative 
expansion of the three- and four-point correlation functions entering the equations 
of motion for <I>(r, i) and ^(r,t) beyond the approximations shown in equations ( pO| ) 
and (^). Their relevance to a microscopic evaluation of relaxation times can be 
understood from the discussion given by Kadanoff and Baym pG|, showing that 
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the analogue of the coUision integral for a normal Fermi or Bose fluid in the Born 
approximation is given by the second-order perturbative expansion of the self-energy. 
In such a viewpoint one may surmise that the three-point term, i.e. and its 

complex conjugate, should give an important contribution to damping of condensate 
motions via correlations with the non-condensate fluctuations. This process should 
exhibit a strong temperature dependence from Bose factors favouring scattering into 
non-condensate states which are being progressively occupied as the temperature 
increases towards its critical value. 
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Appendix A. Comparison with the Hartree-Fock-Bogolubov-Popov 
appr oximat ion 

We show in this Appendix that the RPAE equation ( p^ for the normal modes of 
the condensate reduces to those obtained in the HFB-Popov approximation if we 
choose the excited states of the Gross-Pitaevskii operator for the reference gas of 
the condensate and we neglect the coupling with the density fluctuations of the non- 
condensate. 

The RPAE equation for the condensate can be put in a more explicit form by 
the use of equation for the susceptibility into equation (p^). This leads to an 
expansion of the condensate density change into normal modes: 

<5n,(r,c.) = J2 [K,ir)cb,ir^-\u;)f!j-\u;) - <i>,,(r)^;(r)f7;.+) (c.)/]+)(c^)] (Al) 
where 

r!f^H = l/(a;T(e, -/i)+*7y) (A2) 

with ?7 = 0"*". The coefficients /j ^ and /j^"* in these equations give the selection rules 
for the excitations and are defined as follows: 

ft\^) = I d'r''i>,,{r')cj,*ir') {U^{r',u;) + aSn,{r',u;)) (A3) 

and 

f^+\co) EE / dV'$:^(r')0,(r') {u;ir',co) + a6n,ir\u)) . (A4) 



By substituting equation ( [Al| ) into equation ( |A3|) we finally obtain a set of algebraic 
equations for the coefficients, which represent the RPAE equation expanded on a given 
basis: 

ft^=B^-'>+aY,(ni-^fi-^Ro,oH-ni+^fi+'>R,,oo) , (A5) 



and 
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These involve numerical coefScients defined by 



B 



(-) _ 



EE / dV^*(r')t/p=(r',c.)$(r') 



B^^ EE I d?r' <i>*(r')f/p^(r',c.)0,(r') 



and 



Rabcd 



(A7) 



(A8) 



(A9) 



In equation (|AS| ) we have used the notation <&eq(r) — ?!'o(r)- A similar treatment can 
be given in the general case when the no n-co ndens ate is also present. 

In order to compare with equations (A5) and ( [A6| ), we start from the Bogolubov 
equations in the Popov approximation |l§|, |l3t]: 



(AlO) 



with L'^^ being the Gross-Pitaevskii operator defined in equation (|3^). Extending 
Esry's method at finite temperature , we expand the Bogolubov modes ux and vx 
on the reference gas basis 



where 0/j are defined as the excited states of the Gross-Pitaevskii operator, 



(All) 



(A12) 



By substituting the expansion (All) into the Bogolubov equations ( A10| ) and using 
the orthogonality properties for 0^ (h ^ 0), we find the equations for the coefficients 



f/f)andFi^): 



« E {^oojh Ui^^ - Rohoj V^^^) = (C.A + e, - ^i) T/. 



(A13) 



For Rijki we have chosen the eigenfunctions to be real, consistently with the Popov 
approximation; moreover, uj^^^ and V^^^'' stand for Uhit-^ = lox) and Vh{LO = uj\). 
Equation (A13) agrees with the RPAE equations (A5) and (A6), when we set 
U.ico) = /j") n\-\u;) and V.iio) - /+\u;) n\+\u;). 

In this way we have also checked that in the limit of zero temperature our results 
are compatible with those of Esry [Q, apart from terms of order 1/N where N is the 
total number of particles in the system. 
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